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Abstract. The simulation of particle cascades initiated in the atmosphere by ultra 
high energy cosmic ray particles involves the generation and propagation of a huge 
amount of particles. As it is unpractical to follow every particle to its end, particles 
below a certain energy (Ecut) are discarded from the simulation. In this article we 
study in detail the influence that this cut has on the total energy deposited in the 
atmosphere by the particle cascade in AIRES simulations. The energy deposit is 
directly related to the emission of fluorescence light and is critical for the accurate 
simulation of shower signals in fluorescence detectors. Not correcting for the discarded 
particles introduces a bias on several shower observables related to the energy deposit 
that can range from 3 to 30% or more depending on the Ecut value used. A prescription 
for the correct treatment of these particles is proposed, and the resulting corrections 
to the total energy deposit are addressed, including a new universal parametrization 
of the mean energy deposit per particle. The low energy cut is introduced in the 
simulations to reduce the required CPU time per shower at the expense of simulation 
accuracy. We find that a 0.4 MeV cut for electrons and 0.9 MeV cut for gammas is an 
adequate compromise, and that the proposed prescription is capablable of removing 
the bias introduced by this cut. The prescription is independent of the energy cut 
value and can be used to correct and compare simulations made with different energy 
cuts. 
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1. Introduction 

When a Ultra High Energy Cosmic Ray (UHECR) hits the atmosphere a cascade 
of particles is generated. An important fraction of its energy is deposited in the 
atmosphere as ionization of the air molecules and atoms. A tiny, known fraction of 
the total deposited energy is re emitted during the de-exitation of the ionized molecules 
as fluorescence light that can be detected by ground telescopes. This phenomenon is 
the basis for the cosmic ray fluorescence detection technique pioneered by the Fly's 
Eye experiment |T| and implemented in Hi-Res [2] and the Auger Observatory j3] to 
determine the energy of cosmic rays. 

In this technique, one of the methods to determine the primary energy [1] is to 
convert the light emitted by the shower to the longitudinal profile of charged particles in 
the shower (using the photon yield per particle [5] [6]) and then estimate the calorimetric 
energy integrating the charged particle profile using the mean ionization loss rate, 
obtained from air shower simulations. The estimation of the energy deposited on the 
shower simulation is thus of central importance for the reconstruction and analysis of 
the measured signals. 

The estimation of the aperture and exposure of the detector also relies heavily 
on the correct estimation of the energy deposit. For the simulation of the signal a 
shower would generate on the detectors, the energy deposit is used to estimate the 
photon emission that is then propagated to the detectors taking into account all the 
atmospheric effects. The amount of light arriving to the detector is then used as input 
to simulate the detector response. 

As it is unpractical (if not unfeasible) in UHECR simulations to follow all the 
particles to their rest, a low energy cut is made. Particles with energy below the cut 
value (Ecut) are n °t tracked any more by the simulation program and the fate of the 
energy they carry is not determined. A fraction of this energy might end deposited by 
ionization in the atmosphere and would thus contribute to the energy deposit. This 
introduces a bias on all the observables related to the energy, like the total energy 
deposit, the total electromagnetic energy and the mean energy deposit per particle, 
parameters usually used as input parameters in the simulation and reconstruction of 
UHECR showers signals. The introduced bias has a dependency with Ec u t value, making 
comparison of simulations made with different cuts difficult. 

In this article we present an algorithm to estimate and correct for the bias 
introduced by the discarded low energy particles in AIRES simulations [7], following 
a prescription similar to the one presented in [8] and [9J for CORSIKA simulations. The 
algorithm is tested on a 2 x 10 4 shower library described in section [2j paying special 
attention to the effect this correction has on the shower energy deposit and the mean 
energy loss per particle. To this end, in section [3] we make a careful study of the energy 
balance in the simulation and show the relevance of the discarded low energy particles 
and in section H] we show how to include the contribution of these particles to the shower 
energy deposit. In section [5] we show the effect this has on the average energy deposit 
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per particle, and present a corrected universal parametrization with shower age. 

The low energy cut is introduced in the simulations to reduce the required CPU 
time per shower. Increasing Ec u t value reduces the amount of CPU time per shower, 
at the expense of increasing the total amount of energy discarded from the simulation, 
making the necessary correction more important. We address the influence of Ec u t value 
on the CPU time and on the amount of discarded energy in section [6j 

2. About the simulations 

The quantitative results presented in this work are based on a particular but 
representative set of AIRES simulations. The AIRES code has been extensively used 
by many scientist around the world for the past ten years, and has become one of the 
standard simulation codes in the field. 

AIRES includes the most important processes that may undergo shower particles 
from a probabilistic point of view. For the estimation of the shower energy 
deposit, electrodynamical processes play a central role. Pair production,electron- 
positron annihilation, bremsstrahlung (in electrons, positrons and muons), muonic pair 
production, knock-on electrons (5 rays), Compton effect, photoelectric effect, Landau- 
Pomeranchuk-Migdal (LPM) effect and dielectric suppression effect are all taken into 
account. 

Collision and bremsstrahlung effects have a low energy divergence and can not be 
simulated to zero energy, since the amount of emitted particles would also diverge (not 
to mention the complexity of simulating, for example, electron capture by an atom or air 
molecule). In all simulation codes a threshold energy is imposed, below which emission 
due to these processes is not longer simulated. In AIRES the emission of 5 rays and 
bremsstrahlung radiation is limited to 1 MeV and 0.1 MeV respectively. The particle 
energy loss due to the emission below these thresholds is treated with the continuous 
energy loss approximation during particle propagation, using a parametrization taken 
from GEANT3 simulations. 

These emission energy cuts are independent of and must not be confused with the 
simulation energy cuts studied in this article, that set the energy below which particles 
are discarded from the simulation. In AIRES, different Ec u t values can be set for 
different particle species. The ones more important for the energy deposit study made 
in this work are electron/positron and gammas Ec u t values. 

The shower library used in this work was generated at the in2p3 computing center 
[TO] and has the following characteristics: 

Series name: AMgeLyonExtDvezpShb, generated by Sergio Sciutto 
Hadronic models: QGSJET-II and SIBYLL 

Primary particles: Proton, Iron and Photon (With MAGICS Preshower) 
Log w (Energy): 17.5, 18, 18.5, 19, 19.25, 19.5, 19.75, 20, 20.25, 20.5 
Zenith (deg):0, 18, 25, 32, 36, 41, 45, 49, 53, 57, 60, 63, 66, 70, 72, 75, 78, 81, 84, 87 
Thinning Energy Rel.: 1.0 x 10 -6 Thinning W. Factor: 0.1 
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Gamma Cut Energy = 0.9 MeV 
Electron/Positron Cut Energy = 0.4 MeV 
Muon Cut Energy = 2.5 MeV 
Meson Cut Energy = 4.5 MeV 
Nucleons Cut Energy = 95 MeV 

Thinning refers to a statistical sampling procedure to reduce computing time by simulating 
only a representative sample of all the particles. Thinning and thinning related parameters 
are explained in detail in the AIRES user manual [7]. 

Note that Ec u t for gamma and electron/positron are low enough to determine 
unambiguously the fate of the particles. Gamma below 900KeV will not be able to 
generate pairs, and will deposit all their energy. Electrons generated will also deposit 
all their energy until captured, and positrons will annihilate producing 2 or 3 gammas 
with energy below 1 MeV and thus unable to generate new pairs. 

The results presented in this work exclude photon primaries to make figures clearer 
and keep the discussion simple. The algorithms presented in this article are perfectly 
applicable to photon (or electron) primaries but purely electromagnetic showers have 
qualitative and specially quantitative differences with hadronic showers that make it 
cumbersome to present all the results together. An article focusing in these differences 
is in preparation. 

3. Shower Energy Balance 

During the evolution of a UHECR shower the cosmic ray primary energy is distributed 
among all the particles on the cascade. Part of this energy is deposited in the 
atmosphere, part is lost to the neutrino channel and part arrives at ground level as 
kinetic and rest mass energy of the particles. In simulations, a part is also carried away 
by the particles that fall below the low energy cut and are removed from simulation. 
To check the AIRES simulations and the treatment of the low energy particles in this 
work, we devote this section to the energy balance in the shower. 
In AIRES, the following information is readily available: 

Ez) ep : Total energy deposited on the atmosphere by all particles in the simulation 
(AIRES table 7993). 

Eois'- Total kinetic energy of all particles discarded by the low energy cut (AIRES 
table 7793). 

E v : Total Energy lost to the neutrino emission channel (AIRES .sry file). 

Ecround- Total kinetic energy of all particles reaching ground level (AIRES table 

1793). 

With this information, we can calculate the total energy of the shower as 

SimpleE To tai = E Dep + E Dis + E v + E Ground (1) 
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Figure 1. Histograms showing the contributions to the total energy in the shower for 
our simulation sample. Eu ep :Total energy deposited on the atmosphere by all particles. 
EGround'- Total kinetic energy of all particles reaching ground level. E v : Total 
Energy lost to the neutrino emission channel. _Eoi S :Total kinetic energy of all particles 
discarded by the low energy cut. SimpleETotai is the sum Eoep + Eois + E v + Ecround 
and constitutes the simplest way of computing the total energy on the shower. 



The sum of these four components is the simplest energy balance of the simulation 
that can be made. Figure [I] shows the distribution of the contribution of the different 
terms of[]]in our simulation sample. Epep and Ecround have a strong dependence with 
the shower zenith, since inclined showers traverse a longer distance in the atmosphere 
and deposit more energy, producing the long tails in the distributions. The sum of these 
two components account for 90±5 % of the primary energy. 

Neutrino energy (E u ) accounts for an average 3.4 % of the primary energy in our 
sample, and also have a strong dependence with the shower zenith. In extreme cases it 
can reach up to 10 %. 

Figure [T] also shows that particles discarded by the low energy cut, usually ignored, 
can have an important contribution to the total energy. Even using a low Ec u t value 
like the one used in our sample simulations, they represent an average 5.5 % off the 
total energy. 

The four terms in[T]do not account for all the primary energy but for 99.35 ± 0.4 % 
of it . Some showers are still missing about 1% and there is an overall 0,65 % negative 
bias. This unaccounted energy is carried away by the rest mass of the discarded low 
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Figure 2. Histograms showing the contribution of the rest mass of discarded particles 
and ground particles to the total energy in the shower for our simulation sample. 
Edis rest mass -Total rest mass energy of all the discarded particles. E Groundrestmass : 
Total rest mass energy of all particles reaching ground level. ETotaV- Total Energy, 
calculated as SimpleExotai + Ejji S restmass ~\~ EQ roun d res t mass . 



energy particles, that was not included in the E Dis term of [TJ 

A correction for the rest mass of low energy particles is not straightforward. In 
AIRES not all particle information is saved in the simulation. Low energy particles 
are classified in "gammas" (table 7001), "electrons" (table 7005), "positrons" (table 
7006),"muons +" (table 7007),"muons -" (table 7008), "other charged" (table 7291) 
and "other neutral" (table 7292). 

The particles in "other charged" and "other neutral" are assumed to be pions, 
the lightest unclsassified particle. The contribution to the total energy from this 
miscellaneous groups is of the order of 0.4% when we choose this mass. This is a lower 
limit, as some of these particles could be nuclei, nucleons or kaons. The distribution of 
this correction can be seen in figure [2] labeled as E Disrestmass . 

Not all low energy particles are "created" at the expense of primary energy. Most 
low energy electrons are taken from the atmosphere in Compton or knock on collisions. 
The rest mass of low energy electrons is discarded as a first approximation. 

For completeness, ground particles rest mass will also be included although it was 
found to contribute in average with less than 0,01% to the primary energy. Ground 
particle classification in AIRES is more complete, including classes for protons, neutrons, 
pbar, charged pions and charged kaons (tables 1001 to 1293). Problematic classifications 
are again "other neutral", "other charged" and "nuclei" (all empty in our simulation 
sample). An average rest mass of 1 GeV is suggested for this last three classes, providing 
a lower limit. As in the low energy particles case, ground electrons rest mass is not 
included. The distribution can be seen in figure [2] labeled as Ecround rest mass- 

With the introduction of the rest mass correction the total energy can now be 
written as 

Exotal E[) e p ~\- Ef)i s -\- E v -\- EQ roun d -\- Eui s rest mass EQ roun d rest mass (2) 
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After applying this correction we get a mean total energy of 100.2 % of the primary 
energy, with a standard deviation of 0.2 % (see top right of figure [2]). This results 
confirm that assuming "other particles" to be pions is a good lower limit and that a 
more detailed particle classification is not necessary as this will not change the results 
presented in this article in an appreciable way. The total energy of the shower is expected 
to be always slightly higher than the primary energy, as the medium provides a small 
amount of energy to the shower each time a particle breaks an air nucleus. We can see 
that although considering "other particles" to be pions provides a good lower limit,there 
are still some showers below a 100%. A study on this small "excess energy" might be 
interesting from a theoretical point of view and is planned for a future article. 

4. Energy Deposit off Low Energy Particles 

We have shown in the last section that low energy particles constitute an appreciable 
percentage of the primary energy. Since a fraction of this energy can be deposited in 
the atmosphere, the regular energy deposit estimation in the shower simulation must 
be corrected. 

AIRES provides tables with the longitudinal development of several observables, 
among them the total energy deposit (table 7993) and the energy of the particles 
discarded by the low energy cut {E Disx ), discriminated by particle type (tables 7501 to 
7892). 

Not all the particles lost due to the low energy cut will deposit all their energy, and 
care must be taken when correcting the total energy deposit. Neutrinos from the decay 
of the discarded low energy muons, for example, will not deposit their energy in the 
atmosphere. Discarded low energy positrons, on the other hand, will surely annihilate 
with an electron on the air, giving gammas that will deposit almost all their energy. 
Low energy Pions and other Hadrons might undergo even more complex interactions, 
with an a priori unknown amount of deposited and missing energy. 

An accurate enough correction can be made by adding the contribution from 
electrons, gammas, muons and other particles (assumed to be pions) using the results 
from [8j in which the fraction of "releasable" energy was calculated for each particle 
species using Geant4 simulations: 0.997 for gammas and neutral pions, 0.998 for e+/e-, 
0.425 for muons and 0.46 for charged pions. The prescription suggested to calculate the 
amount energy deposited by the discarded particles {E^ OSit ) is 



The distribution of this correction for our simulation sample and the contribution 
from each term in [3] can be seen in figure [3l It can be seen that the contribution from 
muons and pions is negligible when compared to electrons and gammas. More than 98 
% of the low energy particles are electrons and gammas, that will in the end deposit 
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Figure 3. Histograms showing the contribution of the discarded particles to the 
total energy deposit in the shower for our simulation sample. The top row shows the 
different terms of [3j The bottom row shows the contribution of the kinetic energy and 
the rest mass of the discarded particles and the % correction these represents to the 
total energy deposit of all the charged particles. 



almost all of their energy, making the results insensitive to the details of the contribution 
from other particles. For the Ec u t values used in our simulations, low energy particles 
deposit on average 5.4% of the primary particle energy. 

A second order correction can be made including the contribution of the rest mass 
of discarded unstable particles. Unstable particles that where discarded will eventually 
decay or annihilate, giving electromagnetic particles that will deposit their energy 
producing ionization and more fluorescence light. At this level of detail, the contribution 
of energy from the medium (the air) that is in fact "releasable" energy must also be 
considered. In the annihilation processes, the rest mass of the antiparticle is added to 
the shower energy pool, and part of that energy is also capable of producing ionization. 

Being this a small correction, and following the suggestion in j9] the "depositable" 
energy was assumed to be 1/3 of the rest mass for particles more likely to decay (muons 
and pions) and twice the rest mas for positrons, considered to annihilate and produce 
gammas that will deposit all their energy. As stated previously "other" particles are 
assumed to be pions. The suggested prescription to treat the energy deposit from the 
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rest mass of the discarded particles {E^° r s ^ t mass ) is 

E D?s°rfstrna SS = + ^ n N Dis ^ + 2m e X Dis ^ (4) 

The number of discarded low energy particles (Noi Sx ) is available in AIRES tables 
7001 through 7293. In the simulations made for this work, the correction due to the low 
energy particles rest mass is around 1% of the primary energy, and is almost entirely 
due to the annihilation of low energy positrons (See figure [31 

The total correction due to the low energy particles is obviously dependent on 
Ecut, and this will be addressed in section [6j For the simulations made on this work, 
it has a mean value of about a 7.85 ± 0.33 % of the total energy deposit and is almost 
independent of the shower details. When compared to the energy deposited by the all 
the charged particles above E Cu t (the quantity usually used for shower fluorescence light 
simulations) this correction represents a 8.3 ± 0.3 % increase (figure [3]). 

As the energy deposit is more or less proportional to the emitted light, not 
correcting for the low energy particles introduces a proportional negative bias in the 
signal simulation producing dimmer showers that are harder to detect, reducing the 
detector aperture and that will be reconstructed as events of lower total energy. The 
estimation of the mean energy deposit per particle used in estimation of the primary 
energy in the reconstruction of real detector events is also proportionally affected. 



5. Mean Energy deposit per particle 

The determination of the primary particle energy using the fluorescence technique is 
considered a calorimetric measurement, in which the atmosphere acts as a calorimeter 
where most of the primary particle energy is deposited through ionization. 

The usual procedure [2] [I], consists in converting the measured light profile N ph (X) 
in a region AX to a number of charged particles profile N c h(X) 

Nch[x) = 1 (8) 

AX x Y g geo x g det x g atm 

using the photon yield per electron per g/cm 2 of traversed air Y and the 
corresponding geometrical, detector and atmospheric related correction factors . 

This profile is integrated to get the total number of charged particles and multiplied 
by the mean energy deposit per particle per g/cm 2 < a e ff > 

E em =< a eff > [ N ch (X)dX (6) 



To estimate < a e ff >, the weighted average of the mean energy deposit per particle 
must be calculated using Monte Carlo simulations, this is 

Ground _ , v » 

~ N ch (Xi)AXi chK l> E T ° tal 
e// iaround ~ N T ° tal AX { ) 

i=l 
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Where E dep (Xi) is the energy deposited in a region AXj by the N ch (Xi) charged 
particles crossing the same region. 

Using the prescription for the energy deposit of the discarded low energy particles 
discussed on section HJ the computation of E]^ 1 is straightforward, and to a certain 
extent independent of the simulation Ec u t value as it will be shown in section [6j 

rpTotal 771 1 rpDeposit . rpDeposit /o\ 

^Dep — -^Dep ~r -^Dis 11 Dis rest mass \°) 

The total number of charged particles N^ 1 is not independent of E Cu t, but it is 
possible to apply a correction to render it independent of the low energy cut using the 
parametrization of the charged particles energy spectrum at shower maximum from 

N° ch (Xi) = NS*{Xi)/(l - 0M5E cut /MeV) (9) 

for our simulations, the Ec u t value for electrons is 0.4 MeV and the resulting 
correction factor is 0.982. 

Historically, the standard way of calculating the mean energy deposit per charged 
particle was to divide the total energy deposited by the total number of charged parti- 
cles. For our simulations, the mean energy deposit value obtained in this way is 2.195 
MeV/ g / cm 2 / particle, in excellent agreement with the historical values used Hi-Res and 
other studies [I]. 

If we consider the energy that would have been deposited by the discarded low 
energy particles E^ ostt as defined on[3l the mean value is 2.35 MeV/ g / cm 2 /particle, 
a 7 % increase. Including the rest mass energy correction E^° r s ^ tmass as proposed on 
H]the mean value is 2.375 MeV / g / cm 2 /particle, a 8.2 % increase. This is a correction 
comparable to the correction found in section H] for the total energy deposit, as expected. 

Note that [6] implies that any change in < a e ff > will have a direct impact on 
the estimation of the primary energy. Simulations made with higher Ec u t values are 
expected to scale up this correction accordingly. 

The mean energy deposit is independent (within 0.5%) of the hadronic interaction 
model used and of the primary type, as shown in figure HI This comes to no surprice 
since the energy deposit depends mostly on the details of the low energy electromagnetic 
models. Note that the spread in the distribution of the mean energy deposit for our 
sample is only 2%, showing little dependence with the other shower parameters . 

The mean energy deposit per particle has some dependence with the stage of shower 
development (the shower age, see [Hi), as can be seen in figure [51 This curve and the 
fluorescence yield can be used to convert a Fluorescence photons longitudinal profile to 
a Number of Charged Particles profile, using [5] 

When the shower is young and the number of particles is still relatively low, shower 
to shower fluctuations dominate the dispersion in the energy deposit. Above age 0.7, 
when shower to shower fluctuations are no longer an issue, this curve can be considered 
universal within 1% since the mean energy deposit per particle is virtually independent of 
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<E DBp > [MeV/g.crn 2 /Particle] % diference to Mean Energy Dep % diference to Mean Energy Dep 



Figure 4. Mean Energy Deposit per particle per g/cm 2 of air. Left: Values obtained 
with the full, the simple and without correction for the discarded low energy particles. 
Center: Difference to the mean with by primary particle type. Right: Difference to 
the mean by hadronic model used in the simulation. 



the primary particle, its energy, the shower geometry and the hadronic models considered 
in this article. A parametrization of this curve can be achieved using the function 

£%g*(age) A 
< " e//(age) >= iV* (age) AX = (B + age) c + D + (10) 

age = 3/(l + ^p) (11) 

where A=0.9921, B=0.67, C=9.7878, D=2.1821, E=0.1656 and X max is the depth 
of the shower maximum in gj cm 2 . 

6. Dependency with the Low Energy Cut value 

The algorithms presented in this article to estimate the contribution of the particles 
discarded by the low energy cut are independent of the cut value, and have been tested 
up to 100 MeV. 

The higher the Ec u t value the higher the amount of discarded energy, and the 
impact of the corrections presented in the previous sections on the shower observables. 
To illustrate this, figure [6] (left) shows the dependence of the total amount of discarded 
energy for a set of 10 showers from 1 EeV, 60 deg zenith proton primaries. It can be seen 
that as E Cu t is raised, the amount of discarded energy varies accordingly but the sum 
of the various energy terms of [2] remains constant within shower to shower fluctuations. 

Figure [6] (right) illustrates how the energy deposit on the simulation is lowered as 
more and more energy is lost to the low energy particle cut. But if we consider the 
energy deposit and the corrections proposed by [3] and HI we see that the total energy 
deposit remains independent of Ec u t within shower to shower fluctuations. This provides 
confidence in the proposed algorithms to account for the low energy particles energy, 
and thus to apply them to enable comparisons between simulations made with different 
energy cuts. 
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Figure 5. Mean Energy Deposit per particle per g/cm 2 of air vs Shower Age Values 
obtained with and without correction for the discarded low energy particles. 



For this example 1 EeV proton with the cuts set identical to the ones used in 
our shower library (0.4 MeV for electrons and 0.9 MeV for gammas) we see that the 
discarded particles account for 6.3% of the primary energy. When all energy cuts are 
set at 0.1 MeV (the minimum value AIRES can handle) this goes down to 3.6% of the 
primary energy. Setting the energy cut at 1 MeV it goes to 10.6% , at 10 MeV to 29 % 
and at 100 MeV to 46% of the primary energy (Figure [7]). In this extreme case the low 
energy particles correction to the energy deposit is around a %100 of the normal energy 
deposit. From this, it is clear that extreme care must be taken in the analysis of the 
results of simulations with energy cuts above 1 MeV when the estimation of the energy 
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Figure 6. Components of [2] for the total energy (left) and of [8] for the energy deposit 
(right) of 1 EeV,60 deg zenith proton showers for different values of Ecut- The total 
energy remains constant within shower to shower fluctuations, while deposited energy 
is traded for discarded particles energy as Ecut rises. The total energy deposit is also 
constant when the energy deposit in the simulation is corrected for the discarded low 
energy particles kinetic and rest mass energy deposit 



deposit is important. 

Lowering the low energy cut value to minimize these corrections have a big impact 
on the CPU time. Going from 1 MeV to 0.1 MeV energy cut doubles the required CPU 
time, while lowering it to 10 MeV reduces it by nearly a factor of 3 as shown in figure 
[3 The energy cuts used in our library are considered a good trade off, with nearly half 
the correction found with 1 MeV low energy cut and only 20% more CPU time. 

7. Conclusions 

We have shown the important influence the low energy cut has on the shower energy 
deposit in AIRES simulations. Not correcting for this effect was shown to introduce 
a negative bias on the total energy deposit that can have an important impact on 
shower reconstruction and shower signal simulations. The introduced bias depends on 
the energy cut value used in the simulation, and goes from around 3% of the primary 
energy at 0.1 MeV energy cut to 30% of the primary energy at 10 MeV energy cut on 
average. 
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This Library O.IMeV 1 MeV 10 MeV 100 MeV 

Low Energy Cut 

Figure 7. Influence of Ec u t value on CPU time and % correction to the Energy 
Deposit. This Library has a 0.9 MeV cut for Gammas and 0.4 MeV cut for electrons 
and positrons 

An algorithm to correct for this bias independently of the energy cut value used 
in the simulations was presented and tested for E Cu t between 0.1 MeV and 100 MeV 
successfully, making now possible to compare energy deposit results from simulations 
made with different Ec u t values. Using this algorithm on a large set of simulations we 
computed and provided a new universal parametrization of the mean energy deposit 
per particle with the shower age, that can be used for the reconstruction of the primary 
energy of UHECR shower detected with the fluorescence technique. 

Finally we studied the dependence CPU time has with Ec u t value and found that a 
0.4 MeV cut for electrons and 0.9 MeV for gammas needs only an average 6.3% correction 
due to the discarded low energy particles, giving an adequate trade-off between precision 
in the energy deposit determination and the required CPU time. 
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